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The potential barrier: tunneling 




Figure I: Tunneling through a potential barrier. 

Assume E < V (classically particle is reflected). Outside barrier solutions to the 
SE are 

u(x) = Ae ikx + Be~ ikx for x < -a, (16-1) 

u( x ) = Ce ikx for x > a, (16-2) 

(16-3) 

where we have omitted the term De~ lkx that corresponds to an incident waveform 
the right. Inside the barrier the SE is 

^(x) = +^(Vb - E)u{x) = k 2 u{x) (16-4) 

with k 2 = j£(V — E). As before, k is the decay constant in the classically forbidden 
region is the decay length) that is associated with the "missing" KE necessary 
to surpass the barrier classically, = V — E. Consequently inside the barrier 

u(x) = Ee~ KX + Fe KX , for \x\ < a (16-5) 

As before, we need to match the solution u(x) and its derivative u\x) at the bound- 
aries. 

• At x = —a: 

Ae -ika + Be ika = Ee +Ka + p^na for y ^g_ 6 ^ 

+ikAe~ ika - ikBe ika = +KEe +Ka + KFe~ Ka for u' (16-7) 

• At x = a: 

Ce ika = Ee~ Ka + Fe Ka for u (16-8) 

ikAe lka = -nEe- Ka + KFe Ka for v! (16-9) 



Massachusetts Institute of Technology 



XVI- 1 



8.04 Quantum Physics 



Lecture XVI 



We are interested in the reflection amplitude r — (or the reflection probability 

|r| 2 = |;f | 2 ) and the transmission amplitude t = ^ (or transmission probability 

\t\ 2 = |^| 2 ) from the barrier. Remember that \A\ 2 determines the incident current, 
and is a free parameter. It is useful to divide the equation for u' by the equation for 
u (or alternatively, match = ^ (In u(x)) directly. Then we write 



At x 



-a: 



+ikAe~ ika - ikBe +ika -K,Ee Ka + K,Fe~ Ka 



J^Q—ika _|_ QfAka 



Ee Ka + Fe- Ka 



(16-10) 



At x = a: 



ik = 



-ikCe 



ika 



-nEe~ 



nFe h 



Ce ika 



Ee~ Ka + Fe h 



(16-11) 



(matching of ^(lnu(x)) = ^y^f at boundaries). 
Now we proceed to eliminate E,F (Eq. 16-11): 



ikEe~ Ka + ikFe Ka = — K,Ee~ Ka + KFe Ka (16-12) 

(K + ik)Ee~ Ka = (K-ik)Fe Ka (16-13) 

E = ^^-Fe 2Ka (16-14) 

K + ik v ' 



Substitute into Eq. 16-10: 



-n^Fe 3Ka + KFe~ Ka 
RHS = ^ (16-15) 

K+ik 

—k(k - ik)e +2Ka + k(k + ik)e~ 2Ka 

(« - i£;)e 2Ka + (« + ik)e~ 2Ka 
-K 2 (e 2Ka - e- 2Ka ) + ikK(e 2Ka + e' 2Ka ) 



k 2 sinh(27ra) + ikn cosh(2/ta) 
k cosh(2/ta) — ik sinh(2«;a) 



(16-17) 
(16-18) 
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Consequently, Eq. 16-10 



[+ikAe~ ika - ikBe ika ] [Kcosh(2/ca) - ik sinh(2/ca)] 

(16-19) 

= [Ae~ ika + Be ika ] [- k 2 smh(2na) + ikn cosh(2na)] 

(16-20) 

= Ae~ lka (+ikK,cosh(2K,a) + k 2 sinh(2/ta) + k 2 sinh(2/ta) — ikn cosh(2fta)) 

(16-21) 

= Be lka (+ikn cosh(2fi;a) + A; 2 sinh(2«;a) - k 2 sinh(2«;a) + ikK,cosh(2K,a)) 

(16-22) 

Ae- lka [{k 2 + k 2 ) sinh(2/ta)] = Be ika [2i/U cosh(2/ta) + (k 2 - k 2 ) sinh(2/ta)] 

(16-23) 



r = ~, (16-24) 

A v ' 

= e -2ika (A; 2 + k 2 ) sinh(2/ta) 

2ikKcosh.(2na) + (k 2 — k 2 ) sinh(2/ta) 



reflection amplitude from barrier. 

To calculate the transmission amplitude we use the continuity of u at x = a: 



Ce ika = Ee~ Ka + Fe +Ka (16-26) 
H ~ tk Fe Ka + Fe Ka (16-27) 



k + ik 

-^-Fe Ka (16-28) 
re + ik v 7 
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We find F from the continuity of u at x = —a: 



RHS = Ee Ka + Fe 



k — ik 



Fe 3Ka + Fe 



k + ik 
Fe 



= Fe" 



k + ik k + ik 

2k cosh(2/ta) — 2ik sinh(2/ta) 



k + ik 



RHS = AT ifea + £e ifea 
= Ae~ ika + Ae~ ika 

= Ae-' lka 

= Ae 



(k 2 + k 2 ) sinh(2Ka) 



1 + 



2ikn cosh(2/ta) + (k 2 — k 2 ) sinh(2/ta) 
(k 2 + k 2 ) sinh(2/ta) 



2i/c/tcosh(2fi;a) + (k 2 — k 2 ) sinh(2/ta) 
ika 2ikK cosh(2 ko) + 2k 2 sinh(2/ta) 



2ikn cosh(2/ta) + (A; 2 — k 2 ) sinh(2/ta) 



Then, 

C _ 2k F 
A A n + ik 
2k 



Ka—ika 



Ae 



— 2ika 



[2ikK cosh(2/ta) + 2A; 2 sinh(2/ta) 



A 2/t cosh(2/ta) — 2ik sinh(2/ta) 2ifc«;cosh(2/ta) + (/c 2 — k 2 ) sinh(2fi;a) 

1 



2/te 
C 



-2ika- 



ik 



2ikKcosh(2Ka) + (k 2 — k 2 ) sinh(2/ta) 



-2ika 



2kK 



2kK cosh(2/ta) — ?(A; 2 — /t 2 ) sinh(2«;a) 



(16-29) 
(16-30) 

(16-31) 

(16-32) 
(16-33) 
(16-34) 

(16-35) 

(16-36) 

(16-37) 

(16-38) 

(16-39) 

(16-40) 
(16-41) 



4 



Figure II: Tunneling through the potential barrier. 
Consequently, we have the results for the barrier 
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2m 



r _ B _ c -2ika -i{k 2 +K 2 )smh(2Ka) 

' A 2fcKCOsh(2fca)— i(k 2 — k 2 ) sinh(2Ka) 



^ £ _ C_ _ g— 2ika 

A 2fcrecosh(2Ka)— i(k 2 — k 2 ) sinh(2fca) 



2fcK 



Since the energy and particle velocity are the same on both sides of the barrier, here 
we have Irl 2 + Itl 2 = 1. 




Figure III: The sinh function. 

Let us look at \t\ 2 

\ t \ 2 = ( 2/c/t ) 2 (16-42) 

1 ' (2A;/t) 2 + (P + /t 2 ) 2 sinh 2 (2/ta) V 7 

where we have used cosh 2 (x) = 1 + sinh 2 (x). Since, sinh is a monotonically increasing 
function, and k = ^y/Vo — E, the transmission is monotonically decreasing with 
barrier height Vq. 

In the limit of small transmission, na ^> 1 (barrier width large compared to decay 
length k" 1 ), we have sinh(2/ta) « (|e 2Ka ) 2 = ±e 4Ka and \t\ 2 -> (pq^) 2 e~ 4Ka . In this 
limit the tunneling probability falls off exponentially with barrier thickness (in units 
of decay length k,^ 1 ). 

— > This exponential dependence explains the extremely wide variation in, e.g., 
lifetimes of unstable nuclei (/is to 10 9 years, corresponding to a variation by a factor 
of 10 22 ). 



Massachusetts Institute of Technology 



XVI-5 



8.04 Quantum Physics 



Lecture XVI 




Figure IV: The transmission through the barrier as a function of decay wavevector k. 




Figure V: In the limit of large barrier height or width, the transmission falls off 
exponentially because the wavefunction inside the barrier is dominated by the expo- 
nentially decaying term. 

Potential well: resonance phenomena 

We first consider scattering (E > 0) 

x < -a : Ae lkx + B~ ikx (16-43) 

-a<x<a: Ee +iqx + Fe~ iqx (16-44) 

x > a : Ce ikx (16-45) 




Figure VI: The potential well. 
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h 2 k 2 = E 

2m 



Instead of going through the calculation again, we note that these equations are 
equivalent to those of the potential barrier (for E < Vo) if we replace n — iq. 
Consequently, we obtain 



-2ika (<? 2 - k 2 ) sin(2ga) 



r — i e - llka y -± '- v ± ' (16-46) 

2kq cos(2ga) - i(q 2 + k 2 ) sin(2ga) v 1 

^ _ e ~2ika 2kq 

2kqcos(2qa) — i(q 2 + k 2 ) sin(2ga) 

For the potential well, in contrast to tunneling through the barrier, the reflection and 
transmission oscillate as a function of parameter 2qa, i.e. as a function of number of 
de Broglie wavelengths y inside the well of size a. In particular, for values 

2q n a = nn — > n integer (16-48) 

9» = ^ (16-49) 

K = = - (16-50) 
q n 

the reflection goes to zero because of destructive interference between the waves re- 
flected at —a and +a. This corresponds to the resonance condition for a Fabry-Perot 
resonator in optics, the phenomenon persists in 3D, and for electrons scattering off no- 
ble gas atoms is called a Ramsaner-Townsend resonance. A very similar phenomenon 
has been observed in collision of ultracold atoms, where the effective depth of the 
interatomic potential Vo can be tuned with a magnetic field, there (and in nuclear 
collisions) it is called a Feshbach resonance). 



Bound states in attractive ^-potential 

What happens for negative energies — V < E < in the potential well? 

We expect discrete bound states, at least if potential is sufficiently deep. Particu- 
larly simple mathematically is a limiting case where we shrink the size of the potential, 
simultaneously making it deeper, such that the product of depth and width is con- 
stant. 

Let Vo — > oo, a — > such that a • V — const = A > 0. We then obtain the 
attractive delta potential V(x) = —\5(x). We are interested in bound states: E < 

• Define 

h 2 n 2 

=0-E = -E=\E\, k>0 (16-51) 

2m 1 " v ; 
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Figure VII: If the potential well is sufficiently deep or wide, it can support bound 
states with discrete energies — Vq < E < 0. 



Figure VIII: Attactive delta potential. 



Solutions for x < 0: 



Solutions for x > 0: 



Ae KX + J3^p 

diverges for x — > — oo 



Continuity of wavefunction at x = 0: 



A = D 



(16-52) 



(16-53) 



(16-54) 



Derivative obeys (Lecture XV) 



u'(e) -u'(-e) = -^Am(0) 



(16-55) 



-2k = 



2m 

w 

2m 



A 



(16-56) 
(16-57) 
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m 

«i = — 2 A 



(16-58) 



2,^2 



#1 



ft 2 m 2 



A 



m 



A 2 



(16-59) 



2m "' 2mh i 2h 2 ' 
Binding energy for attractive 5-function. The 5 potential supports 






Figure IX: Comparison of bound states as the potential evolves from a very deep 
to a very shallow potential. In the very deep potential, like in the infinite well, the 
wave function oscillates sinusoidally inside the well, and decays exponentially in the 
forbidden region. In the very shallow potential, the wavefunction is is mostly located 
in the "forbidden" region outside the well. 



exactly one bound state of energy E = —i^t- For a finite-size well, this result 



corresponds to the limiting case of a weak potential that supports only one 
bound state (V <C 



2h z 



with energy E = —^§tVq. 
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Figure X: Solutions in different regions. 

Two attractive (5-potentials 

We could proceed as before, or simplify slightly by making use of the fact that the 
potential is symmetric x — > —x, and therefore we expect solutions of definite parity. 
The even solution in the middle region is 2B cosh(Kx), and A = D, which eliminates 
two parameters. 



• Continuity of u: 



• Derivative: 



2B cosh(/ta) = Ae 



2m 

-K,Ae~ Ka - k2B sinh(/ta) = ^\Ae~ 

hr 



2m 

Y 2 



-X - k ) Ae 



2kB sinh(/ta) 



2m 



2 A — k ) 2_Bcosh(/ta) = 2kB sinh(/ta) 



2ma 

— 2 — A — 1 = tanh(Ka) 



h Ka 



(16-60) 

(16-61) 
(16-62) 

(16-63) 
(16-64) 



There is always exactly one solution of the eigenvalue equation (16-64) for even 
parity. From the figure we see that for the bound state na < 2 ™% x , which is where 
the function A ~ ~ 1 intersects zero. On the other hand, since tanh(x) < 1 
need 
E -- 



h Ka 
-K 1, 



we 



2„2 



h^K 



2m 



or n > pA. Larger k means larger magnitude of binding energy 
A If we compare this to the binding-energy in 



2m 



We have pA < k < 

single (^-potential, Ki = SA we see that the particle is more strongly bound in the 



double-well potential. 
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Figure XI: Graphic solution of the eigenvalue equation 16-64. 



Reason. Given the discontinuity in slope due to the potential, it is possible to choose 
a steeper wavefunction (larger k — > larger binding energy) when the two 5-functions 
are close. Variation of binding energy with well separation a: As we decrease a, the 




Figure XII: Comparison of the wavefunction for two different well spacings. If the wells 
are close, for the same wavefunction discontinuity at each 5 function the wavefunction 
outside the two wells can decay faster (larger k), resulting in larger binding energy 
\E\ = h 2 K 2 /2m. 



FT 

Figure XIII: Graphic comparison of the binding energies for large and small separation 
2a between the binding sites. 

binding energy increases from the value given by k = j£ (binding energy of a single 
well attained at a — > oo) towards the value k = attained as a — > 0. Thus 

the binding energy quadruples, the possibility of the wavefunction in a double-well 
system to change so as to decrease the kinetic (and possibly potential) energy is at 
the origin of chemical bonds in molecules. 




Massachusetts Institute of Technology 



XVI- 11 



